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Abstract 

We consider a model of two-dimensional bosons in a harmonic trap with 
translational and rotational invariant, weak two-particle interaction. We 
present in configuration space a systematical recursive method for construct- 
ing all wave functions with angular momentum L and corresponding energies 
and apply it to L < 6 for all N . The lower and the upper bounds for inter- 
action energy are estimated. We analitically confirm the conjecture of Smith 
et al. that elementary symmetric polynomial is the ground state for repulsive 
delta interaction, for all > L up to L < 6. Additionally, we find that there 
exist vanishing-energy solutions for L > N{N — 1), signalizing the exclusive 
statistics. Finally, we consider briefly the case of attractive power-like po- 
tential r^ , k > —2, and prove that the lowest-energy state is still the one in 
which all angular momentum is absorbed by the center-of-mass motion. 
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Considerable attention has been devoted to the fenomenon of Bose-Einstein (BE) con- 
densation in atomic vapors in the past few years The main reason was the fact that 
such systems might form quantized vortices under rotation, hke in superfluid ^He. Many au- 
thors considered the vortex problem in BE condensates theoretically In the meantime, 
the exsitence of vortex state has been experimentally confirmed |]TO|,|n| . 



Persuing an analogy with the fractional quantum Hall (FQH) effect Wilkin, Gunn, 
and Smith introduced a model of weakly interacting bosons in a harmonic trap [^. In the 
absence of interactions, there is a huge degeneracy corresponding to the number of ways to 
distribute L units of angular momentum among bosons. The degeneracy can be lifted by 
diagonalizing the interaction in the restricted Hilbert space of the lowest-energy states with 
angular momentum L. Accordingly, Bertsch and Papenbrock used numerical computation 
to show that the ground-state energy depends linearly on angular momentum L and that the 
corresponding ground state of the repulsive model for L < is the elementary symmetric 
polinomial of coordinates Zi = Xi + ii/i, relative to the center-of-mass, R= 0- 

Recently, Smith and Wilkin proved analytically that these symmetric polynomials are 
indeed exact eigenstates Later, Huang extended the work by Smith and Wilkin to 



the models interacting via arbitrary rotationally and translationally symmetric potentials. 
The main result was that the associated interaction energy still varies linearly with L. In 
addition, a lower bound for the attractive quartic interaction energy was derived. Finally, 
Bertsch and Papenbrock showed |jl5[ that there is a subspace structure that explains why 



certain eigenstates turn out to be simple analytical functions. They succeeded in construct- 
ing basis states but only up to L = 5. The mentioned subspace structure is present in a 
rather wide class of two-body interactions. However, for general interaction, the problem of 
determining the ground state has not yet been solved analytically. 
Our starting point is the two-dimensional Hamiltonian given by 

^ = E(4v.^ + k)+E^(|r.-r.l). (1) 

i=l i<j 

It describes A^ Bose-Einstein particles in a harmonic trap, weakly interacting via the two- 



particle potential v, possesing translational and rotational symmetries. Since the spacing 
between noninteracting energy levels is greater then the two-body interaction strength, we 
restrict the Hilbert space to the degenerate ground-state manifold of non-interacting 
bosons. The single particle state is of the form z"^ exp(— where z = x + iy and m is 
the angular momentum quantum number. The eigenfunctions of the many-body problem 
are 

N ^ 

tp{zi, Z2,..., zn) = ifizi, Z2,..., zn) n exp(--|zin, (2) 

i=l ^ 

where is a homogeneous polynomial of degree L that is totally symmetric under permu- 
tation of particle indices. Suitable basis functions for such polynomials are given by 

1 ^ / 

where set {Ai, A2, . . . , Xq} denotes any partition of L such that I]?=i \ = L and Ai > A2 > 
• ■ ■ > Ag > for g < |jl6|]. The prime denotes the summ over mutually different indices 
ii,i2, . . . ,iq, while the numbers ni,n2, . . . ,np denote the frequencies of appearences of equal 
X's (i. e., the number of particles carrying the same angular momentum). Note that the 
number of distinct monomial terms z^^z^^ ■ ■ ■ zf^ in Bx is given by A^(A^ — 1) ■ ■ ■ (A^ — Q + 
1) / {ni\n2\ ■ ■ ■ Upl), where ni +n2 + ■ ■ ■ + np = q. For simplicity we omit the exponentials from 
wavefunctions throughtout the paper. Owing to the translational and rotational symmetries 
of the two-particle interaction v, we find, for non-negative integers n and m 

V{ZI + Z2TiZi - Z2)2™P(2;3, 2:4, • • • , Zn) = C2miZl + Z2T{Zx - Z2f'^P{Z'i, Z4, • • • , Zn), (4) 

where P denotes an arbitrary polynomial in ^3, 2:4, . . . , z^v variables, and Cn is given by 

_ dr r2"+it;(r) exp(-rV2) 

~ /o°°c/rr2-+iexp(-rV2) ' ^ ^ 

It is evident that the coefficients c„ represent the interaction energy vir) of the relative 
motion of two bosons in the single particle state r" exp(— r2/2) with angular momentum 
n. We are now in position to estimate the lower and the upper bound of the interaction 



y = T,f<j v{\ri-Tj\). Since Cmm < < c^ax where Cmm(max) = min(max){co, C2, . . . , C2[l]}, 
the total interaction energy hes between 

I 2 JCmin < ^ < I 2 JCmax- (6) 

In order to simphfy the calculations concerning the action of the potential V on polynomials 
Bx, let us define symmetric functions of two variables: 

bij{zi, Z2) = ^{z\zi + ziz^), I > J. (7) 

The action of the potential v on bij is given by: 

[?] 

vhij = Y,a\]hu (8) 

1=0 

where i + j = k + l = n, and the coefficients afj satisfy the summation rule J2iLo otfj = Co- 
The first constraint correspondes to the conservation of total angular momentum for a rota- 
tionally symmetric potential, while the summation rule refiects the presence of translational 
symmetry. By using the simple calculations it can be shown that 

2 _ [f 1 

~ on X/ ^'^P^i!j^n-2p,2pi (9) 
^ p=0 



sh=T.(-YlWV (10) 



where 

r+s=q ' 

The coefficients afj represent the two-body matrix element Vijki (see Ref. of the in- 
teraction potential V . The relation (|^) expresses the coefficients afj in terms of interaction 
energies c„ of a pair of bosons having the angular momentum n. It can be verified that 
in the case of constant potential v = f const = C2n and afj reduces to afj = Applying 
potential V onto B^{zi, Z2, ■ ■ ■ , zn) and using the relation (|), we obtain 

VB'1{zi,Z2,...,zn) = coI^]b'1{zi,Z2,...,zn), (11) 
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i. e., the polynomial -B" (-^i, Z2, . . . , z^) is an eigenvector with eigenvalue equal to the lowest 
(greatest) interaction energy, if cq = Cinin(max)- Furthermore, the action of the potential V 
on the product B^Bx reduces to 

VBlBx = B'lVBx, (12) 

for any n and partition A. Owing to the property (p!^ ) it is convinient to choose the natural 
basis for a given total angular momentum L as a set formed by B^ and products B^^^ Bx 
where 2 < K < L, and A stands for the special partition of K; A = {1,1,...,!} = 

V 

K times 

1^ and all other partitions not containing I's i. e., A = {Ai, A2, . . . , A^}, Ai > A2 > 
■ ■ ■ Afc > 2, where i;j=i Ai = i^' and 2 < < A^. For K > N the partition A = 1^ is 
not present and for K = N + 1, N + 2 the remaining set of polynomials Bx is linearly 
independent. For K > N + 2 the above states become linearly dependent Let us 

call the corresponding subspace spanned by above defined Bx polynomials as IB^- For 
example IB2 = span{i?i2}, JB^ = spanjSis}, JB^ = span{i?i4, S22}, = span{i?i5, ^32}, 
IBq = span{i?i6, i?42, i?33, 5222}; etc. The whole L = 6 space can be described as a summ 
of appropriate subspaces: Bf BflBi © B'fB-i ® Bf]B4^ ® B1B5 © Bq. We note in passing 
that these subspaces are not mutually orthogonal. 

Let us now turn to the solution of the eigenvalue problem. Having in mind the trans- 
lational and rotational invariance of the potential V it is obviuos that VBi^^ Bk = 
B^-^VBk and VBk C span{VBK, BiVBk-i, B'IVBk-2, ■ ■ ■ , Bf}. These give rise to 
the following sequence of the subspace relations: 

VB^-^B2 C span{B^,B^-^B2}, 
VB^-^B^ C span{B^,B^-''B2,B^-^Bs}, 

i (13) 

Consequently, the matrix of potential V in the natural basis possesses the block triangular 
form which means that the original eigenvalue problem has been substantionally simplified, 
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i. e., reduced to the eigenvalue problem in the IBk subspaces, ioi 2 < K < L. Hence, we 
have to calculate VBil, VB22, VB32, VB42, ^-833, VB222, etc., then express the result in 
terms of our natural basis and finally project them onto the IBk subspace. Applying the 
potential V onto the state Bil yields: 



N 



VBil{zi,Z2, . . . ,Zn) = VY^ 

i<j 

uHL-l) 



N -2\ . jN-2\ 



'N-2^ 
. L , 



Co 



«n— ^— + CoL{N -L) + ^{N-L){N~L- 1] 



where the coefficients a are given by relation 

Co + C2 



+ ail n 



B. 



21^-2, 



(14) 



a 



20 



-, ail 



Co — C2 



(15) 



In (14) we have explicitly extracted the {zi,Zj) pair of coordinates to simplify the action of 
the potential V. Binomial coefficients with index t symbolically denote the number of terms 
in remaining N — 2 variables. By using the identity: 



B21L-2 = BiBiL-i — LBiL, 



(16) 



we obtain 



VBil{zi, Z2, . . . ,Zn) 



Co — C2 



{N -L + l)BiBiL-i + 



Hence, the eigenvalue is given by: 



AiL = Co 



cq — C2 
4 



NL. 



BiL. (17) 



(18) 



Having find the eigenvalue we are now in the position to find the corresponding eigenf unction. 
By expanding the unknown eigenvector in terms of basis polynomials Bil and Bi^^B^ with 
unknown coefficients, we easily obtain the system of simple recursive relations. Solving these 
equations we finally get the expansion coefficients in which all dependence on the details of 
interaction V , namely on cq and C2, simply drops out! The corresponding eigenvector is: 



^- = E(-m-"(^) =51.(^1-^,^2-^,...,^^-^), (19) 



N' 



N 



N 



in agreement with Ref. |T^. In the same way we find: 

'N -2 



VB. 



22 



«2l + 2<(iV-2) + co| 



B22 + \a^^B;^ + ^^«2254 + 2«^^52n, (20) 



VB 



32 



+ af^B^i + {N- l)a^^55 + 2a^j522i + 2a^j53ii 



B 



32 



(21) 



42 



«42 + «lo(A^ - 2) + a^[J(iV - 2) + Co 



(N -2^ 



\ 2 



5, 



42 



+ 2a^^533 + a^^^si + {N- l)a'^iB^ + ^^sai + 60^^5222 + 2a^i54ii, (22) 



33 



atl + 2c4^{N-2) + c,\ 



(N -2^ 



1 

— ( 
2 



-B33 + -"33-842 



+ o«33551 + 



A^- 1 



C^SS-^e + "30-^321, 



(23) 



VB. 



222 



3a|^ + 3a^^(A^-3) + Co|^ 



(N -?>^ 



1 

— t 
2 



-B222 + -«22-8321 + 



A^- 2 



"22-842 + 20:20-82211 • 



(24) 



In order to express the above results in terms of the vectors of our natural ba- 
sis we multiply systematically all polynomials B\ by Bi in the following sequence: 
Bl; BiBii,BiB2] BiBia, B1B21, B1B3; BiB^, B1B22, B1B211, BiB^i, B1B4; etc. These 
products can be easily evaluated by simple algebraic manipulations. Here, we give a few 
results of these manipulations: 

Bf = B2 + 2Bii, 
BiBii = B21 + 3-Bi3, 

51^2 = -83 + -821, 
BiBr, =5211+4^14, 

-81-8211 = -8312 + 2i?22i + 3i?2i3, 



(25) 



It is obvious that any polynomial can be expressed in terms of natural basis by method of 
succesive substitution. For example, from the first three relations we obtain = B\ — 

Next, we turn to the final step. Projecting out the states B^^B^i and i?2ii onto S22, 
then S5, ^41, ^312, i?2i3 and ^221 onto ^32, and Bq, B51, B321, 52211,-8313,-8214 and -B412 onto 



we obtain 






P{B^) = 


-5(53i) 


= 2522, 




= 0, 




P{B,) = 


-5(54i) 


= 5532, 


5(53lO 


= 2532, 




P(5213) 


= 0, 




P{B22l) 


= -532, 




pm = 


-5(55i) 


= 9542 + 16533 - 185222, 


P{B32l) 


= —542 — 


2533, 


P{B2Hi) 


= -^5(5313) = i(542 + 2533 - 35222), 


P(5214) 


= 0, 




5(5412) 


= 4542 + 


8533 — 95222, 



(26) 

where P means the projection. Note that vectors 542, 533 and 5222 are linearly independent 
if > 4, while for = 3 we have: 5(542) = —2533 + 35222- Upon substitution of results 
(p6|) into relations ( pO| - p^) we find two more, yet not known, eigenvalues: 

(27) 

O 

'N-2 ~ 



A22 = Co^—{N --) + C2^^ + C4— N>2, L = 4, 



A32 = ^ 

^ 16 



17A^-36 + 8 



+ ^{5N-12) + ^{5N + 2), N>3, L = 5. (28) 
8 lo 



The first one in the case of repulsive S potential (cq 7^ 0, and Cn = 0,Wn > 0) reduces to the 
result already obtain in Ref. |]T^. The two corresponding eigenstates are obtained in the 
analogous way as Ail: 
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A22 = NBi - 3Bl 



(29) 



A32 = NB, + 2B,B2, 



(30) 



where we have used compact notation: Bn = J2f^i{zi — Bi/N)^. Few comments are in order. 
The two above eigenvectors have been aheady found in Ref. |T5[ for a wide class of two-body 
interactions. They do not depend on the details of interaction. But, starting with L = 6 
the wave functions and eigenvalues in the subspace spanned by vectors 3^2, B-^^ and -B222 
become dependent on the details of interaction. Generally, for any L it can be seen from 
relations (|22| - p4D that the action of the potential V ow IBl reduces to VB\ = J^u^xB^ 



where the matrix elements have the structure = y^A/^A^^ + P\fj.N + jXfj.- All details of 
the interaction are incorporated in the constants (3 and 7. For large A^, the eigenvalues of 
this matrix can be easily found by expansion in 1/A^. For example, for the weak 6 repulsive 
potential the corresponding 3x3 (A^ > 4, L = 6) matrix can be written as 



/ 



Co 



V 



_ 51N_ I 15 
2 32 4 

TV _|_ 33 
2 ~'~ 8 

9Af _ 45 
16 16 



9TV 
64 



2 



3 
2 

7TV _ 15 
4 16 



9N I 63 
32 32 



TV I 1 
16 8 



1 
2 

-2A^ + 



8/ 



(31) 



The corresponding eigenenergies can be easily expanded in powers of 1/A^ up to 0{1/N): 



A 



(1) 



A. 



(2) 



A, 



(3) 







Co 




K 2 






Co 








Co 





-2A^ + 



27' 



47 

32' 



34; 
30' 



13; ' 

2955\ 

3536/ 



(32) 



For A^ > 10 the error is of order ~ 1%. Having in mind all obtained eigenvalues, it turns 
out that for L < 6 and arbitrary N > L, the lowest energy is given by A^l. In this way 
we analytically support the conjecture of Smith et. al. [131, up to L = 6. Moreover, the 
first eigenvalue from the relation (P^ ) is in a perfect agreement with the numerical result 
for two excited octupole modes [|T^. For L = A^ it is known that the ground state is a 



vortex ^ = Yl^^i{zi — z^.) 1 around z^, = Bi/N P,|T9[- For L > N, we find that there exists 
lower eigenvalue then that given by ([ISD. For example, for L = 6 and = 4 eigenvalue 



Ag^'' = 1.38 is lower then Ai4 = 2. The same is true for L = 6 and = 3, i. e. Ag^'' = 
and Ai3 = 3/4. For = 3 we find that V^A^iAl^) ~ B2p+3q, and therefore, for a given L 
there exists only one eigenstate having positive eigenvalue VBl = [1 — {—1/2)^^^]Bl JlT 



From the above examples, it seems that for sufficiently large L > N there exists the ground 
state with energy lower than A^iv = cqN{N — 2)/4. Moreover, we find vanishing eigenvalues 
for L > N{N — 1). As 6 potential gives vanishing c^'s for n > 0, from relation (^) the 
corresponding eigenvectors are given by 

N 

^0 = Y[{Zi - ZjfijL-N{N-l), (33) 
i<j 

where ?/'j;,_jv(7v-i) denotes any state in the space of total angular momentum equal to L — 
A^(A^ — 1)- We note that these states saturate the lower bound in inequality (H). It is 
interesting to note that the lowest states (|33D are of anyonic (fermionic) type in the sense 



that they vanish when two particles coincide In other words, particles behave like 

the hard-core bosons. This might signalize some sort of statistical transmutation and of 
exclusion statistics [pl|] . 

We should finally point out that for the attractive, weak two-body potential v{r) = rjr^ 
(for A; < and rj < 0, respectively), the coefficient c„ is given by 

^ ^ r(n + l) ' 



and is always greater then 



co = 2^/V(l + ^)- (35) 



This means that the state is the lowest-energy state for fixed L in which all angular 
momentum is carried by the center-of-mass motion Zc- For k = 4 and Cq = Srj the corre- 
sponding energy equals 8(^^^?7 in agreement with the result of Ref. |T^. The Bf remains 
the lowest-energy eigenstate for the attractive 6 potential (cq = — ?7/27r < c„ = 0) and for 
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all attractive two-particle potentials for which cq = Cmin- This generalization is a simple 
consequence of inequality In all of these examples the presence of the uncondensed 
state is not a surprise because of the attractive nature of the two-body interaction 0. 

In conclusion, we have developed an analytical recursive method in configuration space 
for calculating all eigenstates and corresponding eigenvalues for trapped bosons interacting 
via weak, translationally and rotationally symmetrical pairwise potentials. We applied the 
method to states with L < 6, for all N. We have derived the lower and the upper bounds 
for general interaction energy. For repulsive 5— interaction we have proved the conjecture 
of Smith and Wilkin on the ground state structure up to L < 6, for N > L. Moreover, we 
have demonstrated that there exist eigenstates with vanishing energy, for L > N{N — 1). In 
addition, we have considered the case of general attractive power-like two-particle potential 
r^, k > —2, and found ground state and its energy. 

Note added: After this work was completed a preprint by Kavoulakis et. al. p2[ appeared. 
In the asymptotic limit N oo our eigenvalue Ag^'* from the relation ( p2D is in full agreement 
with their findings in the presence of two octupole excitations. 
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